COINCIDENCE SITE MODULES IN 3-SPACE 



MICHAEL BAAKE, PETER PLEASANTS, AND ULF REHMANN 

Abstract. The coincidence site lattice (CSL) problem and its generalization to Z-modules 
in Euclidean 3-space is revisited, and various results and conjectures are proved in a unified 
way, by using maximal orders in quaternion algebras of class number 1 over real algebraic 
number fields. 



1. Introduction 

A lattice in Euclidean space M'^ is a Z- module of rank d whose M-span is M'^. Two lattices 
r and r' in are called commensurate, denoted by i"' ~ if their intersection F n F' 
has finite index both in F and in F'. This definition also applies to the commensurability 
of Z-modules in M'^ of rank m > d. The commensurator of a lattice F is the subgroup in 
GL((i,M) defined by 

comm(r) := {M £ GL{d,R) \ F ~ MF} . 

It is not difficult to see that comm(i^) ~ GL(d,Q), via the map that takes the standard 
basis of M"^ to a lattice basis. More generally, for a subgroup G C GL(d, M), one defines the 
commensurator of F relative to G in an analogous way, denoted by couiuiq^F). Any com- 
mensurator of this form is a function of the equivalence class defined by ~, i.e., commensurate 
lattices possess the same commensurators. One can extend the definition to the affine group 
and subgroups thereof, but that is outside our scope here. 

Of particular interest in crystallography is the commensurator of a lattice relative to the 
group of isometrics, 0{d, M), or to its rotation subgroup, SO{d, M). This has been investigated 
in great detail in connection with the classification of grain boundaries and twins in materials 
science, compare [HIE]- The connection with the algebraic problem we are interested in here 
emerges via the so-called coincidence site lattices (CSLs), see jSllHE] and references therein. 
These are finite index sublattices of a given lattice F that are intersections of the form F n RF 
with an isometry R from the commensurator. In this context, the notation 

(S)OC(r) := comm(g)o(^j5)(r) 

has proved useful, together with the index function U: OC{F) — > N defined by 

U{R) := [F:{Fn RF)] . 

The integer S{R) is the coincidence index of the isometry R, and the image Z'(OC(-r)) is 
called the elementary (or simple) coincidence spectrum of F. It refers to indices that emerge 
from single intersections only. 

More generally, one considers multiple intersections of the form F D RiF n . . . H ReF, with 
each Ri £ OC{F), and defines S(Ri, . . . , Ri) as the index of this intersection in F. Clearly, 
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the corresponding multiple spectra satisfy 

U{OC{r)"') C i;(OC(r)") for m<n, 

and Utot '■= Um>i '^(OC(-r)"^) is called the total (or complete) coincidence spectrum of F. 
In several important examples, compare [SIHZ]) these spectra are equal or stabilize, in the 
sense that the total spectrum is reached after finitely many lattice intersections. Usually, 
there are more multiple CSLs than simple ones with a given index, though this quantity also 
stabilizes in the examples mentioned. 

Following the discovery of quasicrystals in the early 1980s, there was a clear need for gener- 
alizing these concepts to accommodate aperiodic situations. A very natural approach consists 
in extending the setting from lattices to finitely generated Z-modules that are embedded in 
M!^, such as rings of cyclotomic integers in the plane J19",'3| or similar sets in higher dimensions 
Pl I25L OH] . In general, such sets are dense, but commensurability and the commensurator are 
still well defined, all on the basis of group-subgroup indices, see j2] for details. This is the 
approach we shall use here, too. For a given Z-module F, it leads to the determination of 
its elementary coincidence spectrum and the corresponding coincidence site modules (CSMs), 
classified according to their index in F. In this paper, we focus on single intersections, as 
they permit a rather general solution with algebraic methods based on maximal orders in 
quaternion algebras over real algebraic number fields. Single intersections are also the most 
important ones for the applications in physics, as they refer to the grain boundaries of twins, 
which are more abundant than multiple junctions. 

The paper is written in a self-contained way, in order to also reach readers from the applied 
sciences. For this reason, we have not aimed at maximal generality, and we have tried to give 
simple proofs, rather than refer to general abstract results (though the latter are also men- 
tioned at several places). The full machinery will be developed in a forthcoming publication 
0, together with further examples that are beyond our scope here. 

The material is organized as follows. Section |21 reviews the situation of the body centred 
cubic lattice as a motivating example, hinting at the underlying connection to quaternion 
arithmetic. The connection with quaternion algebra is summarized in Section |21 followed by 
the core of the paper in Sections |3 - IHl Here, we develop the precise connection to reduced 
one-sided ideals in maximal orders and the general solution of the coincidence problem in 
the case that both the quaternion algebra and the base field have class number 1. Sectional 
applies the general findings to three specific cases of interest, which are related to the Hurwitz, 
the icosian and the octahedral rings, followed by some open questions. 

2. Example: The body centred cubic lattice in 

Let be the standard primitive cubic lattice in 3-space, and consider Ibcc := Z'^U(t' + Z^), 
where u = |(1, 1, 1)*. It is well known that 

(S)OC(rbcc) = (S)0(3,Q), 

which is clear from F^cc ~ Z^ and the corresponding statement for Z^, compare [21 Ej. Prom 
now on, we shall restrict ourselves to rotations, i.e., to orientation preserving isometrics. In 
3-space, the orientation reversing isometrics are of the form —R with R a rotation. As — 13 
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commutes with all isometries and is a symmetry of Ibcc (and, in fact, of all lattices and 
Z-modules F in 3-space), it is clear that the consideration of SOC(/bcc) (or of SOC(-r)) 
is sufficient to analyze the corresponding coincidence problem completely. Note also that 
S0(3, Q ) is a dense subgroup of S0(3, M). 

The main interest lies in the classification of the CSLs /bcc ^ -R/bcc that emerge from 
R £ S0(3, Q), grouped according to their index 

U{R) := [Ibcc : (Ibcc n iirbcc)] • 



The answer is well known |13|[2] and can be formulated either explicitly in terms of lattices and 
their bases or in terms of quaternions |21EH]- In particular, the elementary coincidence 
spectrum is the set of odd integers, 

(1) i:(SOC(rbcc)) = 2No + l = {1,3,5,...}. 

An interesting quantity is the number of CSLs of /bcc of index m, denoted by f{m). It turns 
out to be a multiplicative arithmetic function, i.e., f{mn) = f{m)f{n) for m, n coprime, and is 
thus specified by its values at prime powers (other than 1): /(2'") = and f{p^) = {p+l)p^~^ , 
for odd primes p. This results in the Euler factors 



f{p^') 1 + p * 1 — p 



-2s 



for odd primes, and in E2{s) = 1. The corresponding Dirichlet series generating function 
reads 

^-(.) = E 4? = n ^.(^) 

(2) rn=l per 

4 6 8 12 

3^ 5^ 7^ 9^ IP 13^ 15^ 17^ 19^ ■ ' 

where V = {2, 3, 5, 7, 11, . . .} is the set of rational primes and ({s) = Yln>i ^ " Riemann's 
zeta function. As we shall see later, this generating function applies to all three cubic lattices 
(primitive, face centred and body centred), which is why we have called it ^cub- 

This Dirichlet series permits the determination of asymptotic properties of /(m), e.g., by 
means of Delange's theorem (though we do not need the full strength of it), see [H Appen- 
dix] for a formulation adapted to the situation at hand. The result is formulated via the 
corresponding summatory function and reads 
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F{x) := /(w) — x~ = ^r7:rT — , as x — > oo, 



where 1/C(2) = G/vr^ is the residue of <^cub(-5) at its right-most singularity, a simple pole at 
s = 2. This shows that f{m) increases linearly on average and, following ^], see Chapter 
18.2 and the footnote on p. 266 of it, one also says that 6m/7r^ is the average size of the 
arithmetic function f{m). 
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Moreover, as observed in one can express ((2) in terms of the zeta function of the Hurwitz 
ring J of integer quaternions jT^] in the quaternion algebra ]H(Q) as 

, . . _ 1 Cj(^/2) _ (,{s/2) 
^ ^ ^'^'^'^ - 1 + 2-- Ci2s) - C„(V2)' 

where Cj{s) and Cj jl-s) are the Dirichlet series generating functions for the one- and two-sided 
ideals of J, see Section 17.11 for details. This indicates that one should be able to re-derive 
the result using the arithmetic structure of the ring J and to generalize it to other situations. 
That is what we achieve in this article, with special focus on results and conjectures made in 
iniE]- For first results on multiple CSLs of /bco we refer to ISHIEZI- 

3. Quaternion algebra and orthogonal groups 

Let us first recall some basic results about quaternions and introduce some notation, see 
|17[ Sec. 57] for background material. Let be a real algebraic number field of finite degree 
over Q . We consider the corresponding quaternion algebra, 

M{K) = K + iK +iK + kK, 

which is a skew field, otherwise known as a division algebra. A convenient form of the defining 
relations for the generating elements 1 (implicit in the above representation) and i,j,k is the 
original one due to Hamilton |16j . 

i^ = j2 = k2 = ijk = -1, 

together with the requirement that K is central in M{K). We note that K is the exact centre 
of the algebra B.{K). 

A single quaternion x = xo-|-ixi-|-jx2 + kx3 is sometimes also written as x = (xq, xi, X2, X3), 
i.e., as a row vector, while x* denotes its column counterpart. The following result is standard, 
compare |16) and references therein. 



Fact 1. M{K) is canonically equipped with the following structural maps and properties. 

(i) the conjugation: — kxg, which is an anti- automorphism on 
the K-algebra M{K), i.e., one has xy = yx for all x,y £ M{K); 

(ii) the reduced norm and reduced trace: nr, tr: IHI(i^) — > K, defined by 

ni{x) = XX = and tr(x) = x + x = 2xo. 

They satisfy the identity nr(x-l-y) = ni{x)+nr{y)+tr{xy), which makes nr a quadratic 
form on the vector space M(K); 

(iii) the Euclidean scalar product associated to the reduced norm, {x, y) := ^ tr(xy) = 
Ylii^iVij which defines the Euclidean metric on M.{K) = and makes {l,i,j,k} an 
orthonormal basis of M(K) (a so-called 4-bein); 

(iv) an orthogonal decomposition M{K) = E[-|_ © ]HI„ into eigenspaces of the conjugation 
map, where 

H+ := {x e Il{K) \ x = x} = K and 

H_ := {xeM{K) \ x = -x} = iK®iK®kK. 
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The elements of H^, are the scalar quaternions and those of ]HI_ are the pure quater- 
nions. For X G ^{K), we call the corresponding components Re(x) := ^(x + x) the 
real part and Im(x) := ^{x — x) the imaginary part of q. □ 

For X S M.[K), the (reduced) characteristic polynomial (and the minimal polynomial for 
X K) is just X'^ — ti{x)X + nr(x). It is preserved under each algebra automorphism of 
M.{K), and the same is true for conjugation, orthogonality, the metric, and the decomposition 
into scalar and pure part. Therefore, any inner automorphism iq-. x ^ qxq~^ (for x G M.[K) 
and q € M.[K)' := M.{K) \ {0} ) leaves both subspaces IHI+ and ]H[_ invariant and induces an 
isometry on EI_ = K^. Hence, q ^ iq defines a homomorphism ]HI(i^)* — > SO(]HI_) where 
SO(M_) ~ S0(3,i^). 

The map iq fixes the elements of the space K + qK, which, for q ^ K, intersects ]H[_ in the 
line lm(q)K, whose orthocomplement W in EI_ is a plane. Hence, EI_ = lm{q)K © W. For 
w G W, we obtain = 2{q,w) = tr{qid) = qw + wq = —qw + wq. This implies the relation 
qwq~^ = qq~^'W which shows that iq, restricted to EI_, is a rotation with axis Im(g) and 
rotation number qq~^ (whose real part is cos((/9)). 

If g G EI_ \ {0}, then qq~^ = —1, hence the eigenvalues of —iq on ]HI_ are —1 on the 
one-dimensional subspace qK, and +1 on its two-dimensional complement W, i.e., —iq is 
the reflection symmetry of ]HI_. with respect to qK. Recall that any rotation of the three- 
dimensional space ]HI_ is a product of two reflection symmetries —iq, —iq' for q,q' £ E[__, 
see |17| Thm. 43.4, p. 105] for details in the generality needed. Consequently, the rotation 
is induced by iqg', which proves the surjectivity of the homomorphism M.{K)' — > SO(]HI_.) 
defined hy q iq. 

The above argument also shows that, for any q G M{K), one has qx = xq if and only if 
X € K + qK, i.e., g or x is central, or lm{q) and Im(x) are collinear (or parallel). 

The homomorphism q ^ iq can be made more explicit as follows: After identifying K^ = 
\K ®]K ®]^K = ]HI_, the general setting permits rotation matrices in S0(3,-R') to be written 
via Cayley's parametrization (which goes back to Euler for the case K = M., compare jl6j ) as 



(4) R{q) ■= Rg = ——{ 2Kiy + 2Xfi k^-X^ + ^^-u'^ -2k\ + 2^iu 

^^^'^^ \ -2k^ + 2\v 2KX + 2fiu K^-X^-fi^+u^^ 

where q = k + iA + j/i + kz^ G M.{K)' , and where K^ is considered to be represented by the 
column vectors (x, y, zY with x,y,z £ K. 

This explicit parametrization is useful in practice. The rotation axis (when Rq is not the 
identity) is given by (A,/i,i/)*, while the corresponding rotation angle if, using that the trace 
of R{q) is 1 + 2 cos{ip) by Euler's theorem, is determined by 

cos(^) = Re(,/g) = ^ = ^2 + ;,2+^2 + ,2 • 

Fact 2. If K is a real algebraic number field, every rotation matrix M G 80(3,7^^) is of the 
form M = Rq for some q G Il{K)' , i.e., the mapping M.{Ky — > S0(3, JC) defined by the 
Cayley parametrization q^ Rq from (0} is onto. 
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Proof. This follows from the above arguments. Alternatively, one can use the parametrization 
(0} for S0(3) = S0(3,M), where it is well-known since Euler and emerges from the canonical 
2 to 1 map from SU(2) — > S0(3), combined with the observation that R{q) = R{aq) for 
arbitrary non-zero a G M. One can then show that a parametrization of the subgroup S0(3, K) 
can be achieved with the quaternions restricted to M{K)'. □ 

Looking at this matrix setting in 4-space, noting that = M{K) = K (B iK (B jK © kK 
and recalling our convention to write quaternions as row vectors, one has the relation 




for arbitrary a £ M{K) and q G M{K)', which reflects the orthogonal decomposition into real 
and imaginary parts mentioned above. In particular, one also has 

(5) lm{qaq^^) = Rqlm{a), 

which will be crucial for our further development. 

4. Quaternion arithmetic 

Consider the real algebraic number field K and let o be the ring of integers in K, compare 
jni for general background material on algebraic number fields. 

Assumption. We shall assume throughout this article that K has class number 1. 

The important consequence of this, for our purposes, is that o is a principal ideal domain 
(PID) and hence that every o-module we encounter will be free and thus has a basis. Also, 
our results permit a simpler formulation this way, though many can be generalized to avoid 
this assumption. 

Next, we need to know what 'integral' and 'integer' mean in M.{K). A quaternion q is 
called integral (over o), when both tr(g) and nr(g) =: |gp are in o, i.e., when it is a root of a 
monic polynomial with coefficients in o. 

The notation |qp instead of nr(q) is quite common and highlights the fact that the reduced 
norm coincides with the squared Euclidean norm of the quaternion under the standard iden- 
tification of the generating elements with the Euclidean 4-bein. However, it has to be stressed 
here that, in spite of the notation, \q\'^ is in general not a square in K. 

Clearly, if q is integral, then so is q. Note that, because of the invariance of conjugation, this 
definition is invariant under all automorphisms o{'El{K). However, due to non-commutativity, 
the set of integral quaternions fails to be a ring. This is the origin of various differences from 
the arithmetic of commutative number fields. 

Following |2S|, we call a finitely generated o-module whose X-span is the whole of M.{K) 
an ideal. An order is an ideal that is also a ring containing 1, see [101 Sec. 8] and j23j Chap. 1, 
Sec. 4] for details. It is called maximal when it is not contained in any larger order of M.{K). 
Note that some authors restrict the term 'ideal' to o-modules whose left and right orders are 
maximal, because their properties are then more similar to those of ideals in algebraic number 
fields. Another term in use is that of a (full) o-lattice, compare |2fl| I23j . 

The following basic properties of an order L are important throughout the article. 
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Fact 3. All elements of an arbitrary order L C El(i^) are integral. Moreover, one has L = L 
and Ln K = o. 

Proof. The first assertion follows from [23J Lemme 1.4.1 and Prop. 1.4.2] or |2U[ Thms. 8.6, 
9.3 and 1.14]. If x G L, one has x = tr(x) — x, so x £ L since tr(x) G o, and L = L. 

Since 1 £ L and L is an o-module, o C L n -fC is clear. Conversely, if x G Lf] K, one has 
tr(3;) = 2x G O and nr(x) = G O, which implies that also x G C □ 

An element a of an order L is an L-unit when G L, and it is a consequence of Fact 01 
that a G L is an L-unit if and only if nr(a) is a unit of o. When it is clear which order is 
meant we shall omit the prefix L and simply speak of a unit. A left ideal of L is an ideal / 
with al C I for every a £ L and it is principal if it has the form Lx for some x G M{K). 
Right ideals and principal right ideals of L are defined similarly. A two-sided ideal of L is an 
ideal that is simultaneously a left ideal and a right ideal of L. 



is an order, as is qCq~^ for any q G M.{K)'; but it is not maximal because + i + j + k) is 
integral and appending it to C generates a larger order. For the applications, we are mainly 
interested in maximal orders that contain C In the literature, such maximal orders are 
sometimes simply referred to as rings of integers. 

Fact 4. If q is an element of a maximal order O that contains C, then 2q G C. 

Proof. The o-basis of C is also a K-hasis of M(K), so q = a + ib + jc + kd with a, b,c,d£ K. 
Since (7, \q, jg and kg must all be integral (as i, j and k are units of O), their reduced traces 
must be in o. They are 2a, —26, —2c and —2d, respectively, so that 2q £ C. □ 

Maximal orders O are certainly not unique, as the application of inner automorphisms 
results in different maximal orders. They may not even be unique up to inner automorphism. 
Orders related by an inner automorphism are said to be of the same type. The class number 
of a maximal order O is the number of equivalence classes of left ideals of O under the 
equivalence relation of right multiplication by elements of IH(i^). Since the conjugation anti- 
automorphism takes left ideals to right ideals, right multiplication to left multiplication and 
O to itself, interchanging the roles of 'left' and 'right' in this definition results in the same 
value of the class number. 

Fact 5. All maximal orders of M.{K) have the same class number, in the sense just defined, 
which is known as the class number of M{K). Moreover, if the class number is 1, all maximal 
orders are mutual images of one another under inner automorphims. 

Proof. This follows from Lemme 1.4.10 and Exercise 1.4.1 of I— ' 

When E[(Er) has class number 1, every left ideal of a maximal order O has the form Oq for 
some q £ O' := O \ {0} (and every right ideal has the form qO). In the case of EI(Q), there 
is only one type of maximal order, and a unique maximal order O that contains £, see |15| : 



Clearly, 



(6) 



C := O + io + jo + ko 



(7) 



O = S = <l,i,j 



i(l + i+j + k)> 
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This is the so-cahed Hurwitz ring of integers ^Sj. It is an index 2 extension of £ ~ Z^. In 
fact, seen as a lattice in 4-space, it is the weight lattice Dl, the dual of the root lattice D^, 
the latter in standard representation with roots of squared length 2, compare jllj . 

Later on, we shall need some information about the relation between O and qOq~^, for 
q & O. It is immediate that 

(8) O = qOq-^ ^ qO = Oq ^ qO = Oq, 

the last characterization following from O = O. When this occurs, qO = Oq is said to be a 
two-sided ideal of O, and q is said to generate a two-sided ideal. All elements of K commute 
with O, and thus generate two-sided ideals. Another example we shall meet below is 1 + i, 
which generates a two-sided ideal in J, but not in our other explicit examples of Section [7| 

Finally, we observe the following counterpart of the relation Lr\K = o from Fact |31 where 
we use the set-valued extension of Re, i.e., Re(S') := {Re(x) | x G 5}. 

Fact 6. If IHI(i^) has class number 1 and if O is a maximal order in M.{K), O has real part 
Re(C') = \o. 

Proof. Assume first that O is a maximal order that contains £ of © as well as the integral 
quaternion ^(1 + i-|- j +k). It is clear that such a maximal order exists, and that C Re(C') 
in this situation. On the other hand, 2q ^ C for all g G O by Fact in this case, and we must 
also have Re(C') C io, hence Re(C') = \o. By Corollaire 1.4.11 and Exercise 1.4.1 of [13, ah 
other maximal orders of M.{K) are images of this one under inner automorphisms, and, since 
the real part is invariant under inner automorphisms, the result follows. □ 

Below, we shall often meet intersections of two (not necessarily distinct) maximal orders. 
These are also known as Eichler orders, see |231 121| for details. 

In this and the following two sections, we develop some general arguments without using 
specific properties of the maximal orders we have in mind for the applications. To this end, 
let K as above be a real algebraic number field with class number 1 whose ring of integers o is 
therefore a PID. The class number property implies that every o-module in EI (i^T) is free (and 
thus has a basis) ^ Ch. 3, Thm. 6.6], since such modules are torsion-free. In what follows, 
L will be an arbitrary order in M.{K) and O a maximal order. Later on, we shall need the 
assumption that O has class number 1, i.e., that all ideals of O are principal. This condition 
will be formulated as a property of EI(Er) and always be mentioned explicitly. 

Let us return to the parametrization of the rotations. It is clear that the choice of (7 in Q 
is not unique. In fact, we can always arrange q £ O' (by multiplying q hy a suitable element 
of o). An element q £ O is called O -primitive if it is not divisible by any non-unit of O. For 
an arbitrary q £ O', one can define 

(9) cont0(g') := lcm{a G O \ q/a G O} , 

which is called the O-content of q, see jl] for details on this concept in a more general setting. 
It is well defined up to units of o due to our assumption that o is a PID. So, q is O-primitive 
when contQ{q) is a unit in o. Alternatively, one can consider the content to be the principal 
ideal generated by cont(p(g), which takes care of the units automatically. 
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More generally, when o is a PID as here, primitivity relative to any order L of M.{K) can 
be defined in an analogous way: q £ L is L-primitive if it cannot be factorized as q = ar, with 
r £ L and a a non-unit in c As with units, we shall omit the prefix L when it is clear which 
order is meant. Also the L-content, cont^(g), for any q G L' can be defined analogously to 
the C-content, and is again unique up to units of o. In particular, primitivity relative to the 
order £ of © is another useful concept. In IS , an element q £ C that is £-primitive in this 
sense is often simply called primitive, and this property is equivalent to the condition that 
gcd(g) := gcd(gQ, g^, 93) = 1; where the gcd is defined as usual, and determined up to a 
unit of o. More generally, one has cont£(g) = gcd(g). 

Proposition 1. Let L be an arbitrary order of M.{K) and let R : q ^ R{q) Cayley's 
parametrization of Then, one has R{L') = S0{3,K). Moreover, any rotation matrix 

M £ S0{3,K) can be parametrized as M = R{q) with an L-primitive element q £ L. This 
gives a multiple cover of S0(3,i^), where the remaining multiplicity is due to R{£q) = R{q) 
with e a unit in o. 

Proof. Since R{Wi{K)') = S0(3, K) by FactH it is clear that, given a matrix M £ S0(3, K), 
q can be chosen to satisfy R{q) = M and to be an L-primitive quaternion at the same time, 
the latter condition due to the freedom to multiply g by a suitable non-zero element of K, 
which commutes with q and does not change the matrix M. 

Clearly, R{eq) = R{q) for any unit e £ o, and eq is L-primitive when q is. In view of 
Eq. ©, and the fact that Re{qaq^^) = Re(a), two quaternions p,q parametrize the same 
matrix M if and only if p~^q is a central element of M.{K), so q = ap with a £ K. When p 
and q are both L-primitive, such a relation can only hold when a is a unit in c □ 

Remark 1 . The relevance of Proposition ^ originates in the observation that 



since Im(£) is the o-span of the standard Euclidean 3-bein in M^. Moreover, S0(3, K) is the 
SOC-group of Im(L) for any order L of M{K), because C and L (and hence also Im(£) and 
Im(L)) are commensurate. 



One of our main goals is to compute indices of CSMs. As we shall see below, this will 
require working out indices of the form [Im(L) : Im(gL)] for elements q £ L'. To achieve this 
in a simple and systematic fashion, it is advantageous along the way, and essential for our 
later proofs in this concrete approach, to work with indices which take values in o rather than 
in Z, and we define these now. In fact, we define these more general indices to be principal 
ideals in c In this setting, the coset-counting indices of interest can later be derived as the 
absolute norms of the corresponding o-indices, see Fact |H1 below. 

For a vector space V over K and a i^-linear map (j) ofV into itself, we use det ((/>) to denote 
the usual determinant of (p. It is an element of K and can be calculated using a i^-basis of V. 
Similarly, detQ (cp) denotes the determinant of cp when V is regarded as a vector space over 



(10) 




5. Index calculations 
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Q and (p as a Q -linear map. It is a rational number and can be calculated using a Q -basis 
of V. In both cases, the result is independent of the particular basis chosen. 
If J C / are full o-modules in V, we define the K-index of J in I as 

(11) [I:J]k := detMc 

where (/> is a linear map that takes an o-basis of / to an o-basis of J and detj^{(j))o is the 
principal ideal of o generated by det;^(i;^). The determinant detj^{(j)) is in o, since J ^ I, and 
changing the bases of / and J multiplies it by a unit of o, so the ideal [/ : J]^ is independent 
of the chosen bases. By considering the inverse of 0, it is clear that I = J if and only if 
[/ : J]^ = o. The X-index is multiplicative: if Ii ^ I2 ^ I3 are full o-modules in V, one has 

(12) [h : h]K = [h ■■ h]K [h : h]K- 

We also write [I : J] = |detQ (0)|, a positive integer, which is the index of J in / in the usual 
sense. 

Lemma 1. For any order L in M{K) and any q £ L' , one has [L : qL]j^ = [L : Lq]j^ = \q\^0. 

Proof. It is sufficient to prove the claim for right ideals qL. The linear map (pq: x ^ qx (the 
left regular representation of q on M(K)) takes an o-basis of L to an o-basis of qL, and a 
straightforward calculation using the basis {l,i,j,k} oiIl{K) gives detj^{4)g) = \q\^. □ 

Our next result is an index formula that provides the link to our problem in 3-space by 
giving the indices of images under the mapping x t-^ Im(x) and paves the way to 3-dimensional 
analogues of Lemma ^ 

Proposition 2. If J I are ideals in M{K), one has the index relations 
[I:J]k = [/ n M+ : J n [Im(/) : Im( J)]^ and 

[I : J]k = [Re(I) : Re( J)]^ [/ n M_ : J n M.]^ , 

where all five indices are principal ideals of O . 

Proof. These identities are both instances of the fact that the density of a lattice (being the 
reciprocal of the absolute value of its determinant) is the density of an orthogonal projection 
of it multiplied by the density of the kernel of the projection. 

For the first identity, choose generators and /Jg of / H EI-|_ and J H and extend them 
to o-bases 

{ao,"i + '"i,"2 + '"2>«3 + ^3} and {^q, (3^ + u^, f)^ + u^, (3^ + u^] 

of / and J; these operations being possible by Thm. 6.16], since ckq and /Jg are primi- 
tive elements of / and J in the sense of p . By construction, {Im(ii]^), Im(n2), Im(u3)} and 
{Im(u;^), Im(?;2), Im(u3)} are then O-bases of Im(/) and Im(J). Note that, due to our con- 
vention, each Im(uj) is the transpose of the row vector Ui without its first coordinate, and 
similarly for the row vectors Vj. Since /nBI_|_ and JnBI_|_ are 1-dimensional, [/nIH[_|_ : Jn]H[_|_]^ 
is trivially evaluated as (/3g/ao)c- If 

X = fiQUo + /xi(ai + ui) + H2{(^2 + U2) + 1^3(0^3 + U3) > 
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one has 

lm{x) = lm{fi^ui + H2U2 + /^3%)> 
so the 4x4 matrix of the hnear map taking the basis of / to the basis of J has the form 



M 



M3 



M 



where M3 is the 3x3 matrix of the hnear map taking the basis of Im(I) to the basis of Im( J). 
Taking determinants of both sides gives 

[/ : J]k = det(M)o = (/3o/ao) det(M3) o = [/ n M+ : J n 11+];,. [Im(/) : Im(J)]^. 

For the second identity, let a and P be generators of the o-ideals Re(/) and Re(J), choose 
Uq and Vq in IH_ with o + Uq £ I and P + Vq £ J, and choose bases {ui,U2,U3} and {^'1,^2,^3} 
of / n EI_ and J n EI_. Then, {a + uo,ui,U2,U3} and {P + vo,vi,V2,V3} are bases of / and 
J, and the 4x4 matrix of the hnear map taking the former to the latter is 

^P/a m 
D Ml 

where M3 is the 3x3 matrix of the linear mapping that takes {Im(u^), Im(u2); 1111(^3)} to 
{Im(W]^), Im(f2), 1111(^3)}. Hence 

[/ : J]k = det(M)o = (p/a) detiMs) o = [Re(/) : Re(J)];^ [/nH_ : JnU.]^, 

which completes the argument. □ 

Since EI4. ~ in our setting, we shall usually write I DK instead of / H H-i- from now on. 

Fact 7. If L' C L are orders in M{K), one has [L : L']j^ = [Im(L) : Im(L')]^. 

Proof. This follows from Proposition |21 since L n K = L' n K = o, hy Fact 01 □ 

Lemma 2. If L is an order of M{K) and q £ L is L-primitive, one has 

qLDK = LqDK = \q\^o. 

Proof It is sufficient to prove qL Ci K = \q\'^o, as the other claim follows by conjugation. 
Observe that the inclusions \q\'^o C qL n K and qL Ci K <Z o are clear. It remains to show 
that a £ qL n K implies |gp | a. We have qa = aq = \q\'^r, where also q is L-primitive, i.e., 
cont^(g) is a unit in o. Consequently, we have {ql"^ \ cont^(ag) = a. □ 

Lemma 3. If L is an order of M{K) and q £ L is L-primitive, one has the index formula 

[Im(L) : lm{qL)]K = [lm{L) : Im(L^)]^ = \q\'^o. 

Proof. This follows from Proposition El and Lemmas ^ and |21 since [o : \q\'^o]j^ = \q\'^0. □ 

In the three lemmas derived so far, and in various places below, there are obvious additional 
identities, due to the freedom to replace q by q. 

To make the transition from if-indices to the usual indices of modules, we need to introduce 
norms from K to Let N := IN^/q] be the absolute norm from ivT to Q , which is given 
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by N(a) = | n^('^)| ^ ' where a runs through the homomorphisms of K into C. It is a 
non-negative rational number for all a G if, and a non-negative integer for a G c 

The absolute norm can also be defined for principal ideals of o by N(aC') = N(q;), and 
is independent of the generator a, since the ratio of any two generators is a unit and has 
absolute norm 1. 

Fact 8. With N as above and V, (j) as in the paragraph preceding Eq. ()11() . we have the identity 
\detQ{(j))\ = N(det^((/))). In particular, [I : J] = N([/ : J]j^) for any ideals I, J of Il{K) 
with J Q I. 

Proof. This is a standard result that can be proved by embedding y in V^© ~ M"*^, where 
Vr is the real vector space V ®k ^ (the realification of V) , Vc is the complex vector space 
V ®K C (the complexification oiV), r is the number of homomorphisms of K into M, s is the 
number of complex conjugate pairs of homomorphisms of K into C, d = r -|- 2s is the degree 
of K over Q and n is the dimension of V. Then, (j) can be extended to an M-linear map of M"'^ 
into itself and the determinant of the extended linear map is the norm of the determinant of 
(j). The 1-dimensional case of this construction is given in |2l Ch. II. 3.1] and the general case 
in dHl Sec. 4]. □ 

This gives the following useful consequences. 

Proposition 3. Let L' C L be arbitrary orders in M{K) and q €z L' . Then, one has 

(i) [L:L'] = [Im(L):Im(L')]; 

(ii) [L:qL] = [L:Lq] =N{\q\^); 

(iii) [Im(L) : lm{qL)] = [Im(L) : Im(Lg)] = N(|gp), provided that q is L-primitive. 

Proof. Part(i) is immediate from Facts [7| and |H1 Part (ii) results from combining Fact |H1 with 
Lemma ^ and Part (iii) from combining it with Lemma |21 □ 

Remark 2. Part (ii) of Proposition IHl clearly displays the correct scaling behaviour of the 
index when q is replaced by aq with a £ o. Taking the scaling behaviour into account in 
Part (iii), it is almost immediate that its extension to general q £ L' reads 

(13) [lm{L) : lm{qL)] = [Im(L) : Im(Lg)] = N(conti(g)) N(|g|2) . 

It is possible [211 to derive generalizations of (|13|) to central simple algebras of dimension 
over fields with a Dedekind ring o, with q i— > lm{q) replaced by the mapping q ^ q — ;^tr(g). 
Our formula would then result from taking n = 2 and considering only principal ideals. 

Let us now look at the coincidence problem for the rank 3 o-module P := Im(L), where we 
want to classify the intersections PCiRP with R £ SOC(-r) according to their indices. From 
Proposition we already know that, for finite index, we may restrict to R = Rq with q £ L', 
and with the identity Q we have 

(14) PnRqP = Im(L)ni?gIm(L) = lm{L) nlm{qLq-^). 
Lemma 4. If L is an order of M(K) and q £ L' , one has the relation 

Im(L) nIm(gLg"^) = lm{L r\ qLq'^) . 
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Proof. The inclusion 5 is immediate, so we have to prove the converse. Consider an element 
n* G Im(L) n Im(gLq^^), where we write for the transpose of u without its first coordinate. 
As lm{qLq^^) = Rglm{L), there is some G Im(L) with u* = RqV^, and we can choose 
X = (a, u) and y = {(3, v) in L, with a, P £ K. Assume for a moment that a — /? G c Then, 
z = a — /3 + y = (a, v) £ L with Im(2;) = Im(y). Moreover, one has 

qzq^^ = q{a,v)q~^ = {a,u) = x £ L 

so that qzq~^ £ Ln qLq~^ with Im(gzg^-'^) = u*. 

It remains to show that a — P £ o. Observe that, since x and y are integral, 

tr(x) = 2a := 7 G o, 
tr(y) = 2/3 := (5 G o, and 
nr(x) - nr(y) = - = (7^ - 5^)/4 G O, 

the last line because nr(y) = ni{qyq~^) = nr(/3, u). We now have 4|(7^ — 5^) with 7, 5 G O, and 
the proof will be complete if we can deduce from this that 2 1(7 — 5). Let vr be a prime factor 
of 2 and suppose that vr^ || 2. Now, Tr*^ f (7 — 6) would imply vr^ f (7 + 6), since j + 6 = ^ — 6 
(mod 2), giving vr^*^ f (7^ — (5^) and contradicting the fact that 7^ — 6'^ is divisible by 4. Hence 
every prime factor of 2 divides 7 — 5 to at least as high a power as it divides 2 and we have 
established that 2|(7 — 5) and thus a — {3 £ o. □ 

The set of L-primitive elements of L is still too big for our purposes. Recall that any 
q £ L' with the property that qL is a two-sided ideal also satisfies qLq~^ = L, and is thus not 
of interest for the coincidence problem (at least as far as counting the CSMs is concerned). 
In fact, the elements of L that generate two-sided ideals give, via the parametrization R, a 
multiple cover of the rotation symmetry group of the module Im(L), the multiplicity being 
due to units of c Consider a general element x £ L' and assume that it is of the form x = aqb 
with a,b,q £ L' and qL = Lq, so that q generates a 2-sided ideal. Then, there are elements 
a', b' £ L' such that x = qa'b = ab'q, and q is then both a left and a right divisor of x. We 
can thus simply call q a divisor of x in this case. Now, x is called L-reduced if no divisor q of 
X other than the units of L satisfies qL = Lq. Again, we shall omit the prefix L when it is 
clear from the context. Recalling Eq. (0), and thus the identity 

Rglm{L) = Im(gLg-^), 

the following consequence of Proposition ^ and Lemma 0] is immediate. 

Corollary 1. Let F = Im(L) and consider a CSM of F , obtained from a rotation matrix 
M £ SOC{F) = S0(3, i^). Fhen, among the SOC{F) -matrices that result in the same CSM, 
there is at least one of the form R = R{q) with q an L-reduced element of L. □ 

Consider now a maximal order O. Our next aim is to relate the conjugate orders qOq~^, 
which are all of the same type as O, to ideals qO of O. For this, we shall need to know 
something about the irreducible factors of elements of O, and the following fact is sufficient 
for our purpose. 
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Fact 9. Let M{K) have class number 1, and let O be a maximal order in M{K). If q ^ O 
is primitive and -k ^ o is a prime factor of \q\^ , then q can be factorized as q = rp, with 
r,p £ O and |pp an associate of n. 

Proof. We translate the proof of |12| Ch. 5, Thm. 2] to this more general setting. 

Since M.{K) has class number 1, the left ideal Oq + On has the form Op for some p £ O 
(the greatest common right divisor of q and vr, determined up to multiplication by a unit). So 
q = rp and n = sp, for some r,s £ O, and vr^ = nr(7r) = Since vr is prime, \p\'^ must 

be an associate of 1, vr or vr^. However, p is not a unit because p = tq + nvr, for some t,u £ O, 
so nr(p) = nr(t) nr(g) + nr(M)vr^ + Trtic{tqu) is divisible by vr. Also, if \p\'^ were an associate of 
vr^, then s would be a unit, making vr an associate of p and therefore a divisor of q, which is 
impossible since q is primitive. Consequently, |pp is an associate of vr. □ 

We now embark on finding a correspondence between qO and lm{qOq~^) when q is reduced. 

Lemma 5. Let Il{K) have class number 1 and let O be a maximal order. If q,r £ O are 

primitive, one has 

Im(gO) C Im(rC') =^ qO C rO and qO ^ rO =^ On qOq~'^ C O n rOr~'^. 

Proof. Since M{K) has class number 1, the right ideal qO + rO has the form dO for some 
d £ O (the greatest common left divisor of q and r). When lm.[qO) C Im(rC'), we have 

lm{dO) = lm.{qO + rO) = Im(5C') + Im(rC') = Im(rC'), 

so, by Lemma 01 and recalling = nr(g), 

m:{d) o = [Im(C') : Im((iC')]^ = [Im(C') : Im(rC')]^ = nr(r) o, 

since d (being a left divisor of r) is primitive. Hence nr(r/d) is a unit in o, so r/d is a unit 
in O and rO = dO 3 qO, establishing the first implication. 

Now assume that qO C rO. Then q = rp for some p £ O. Given any x £ O H qOq^^, put 
y = r^^xr and define 

Iden(y) := {a £ O \ ay £ O} 

(the left denominator of y). This is a left ideal of O and has the form Od for some d £ O. 
Clearly, lden(?/) contains r, but it also contains p, as y = pzp~^ for some z £ O. Consequently, 
r = sd and p = id with s,t £ O, and hence q = \d\'^si. Since q is primitive, |dp is a unit of o 
and hence d is a unit of O. Thus Iden(y) = O, so y £ O and x = ryr~^ £ O Ci rOr^^. This 
establishes the second implication. □ 

Lemma 6. // M.{K) has class number 1, O is a maximal order, and q £ O is reduced, then 
Im(gC>) = Im(CngCg-^). 

We note that, except when g is a unit, it is only the projections into BI_ that are equal: 
qO and O H qOq~^ themselves are distinct (because 1 £ qO only for q a unit). 

Proof. We use induction on the absolute norm N(|gp). The result is certainly true when 
l^(kP) = 1) since then is a unit in o so g is a unit in O and qO = O = Or\qOq~^ . When 
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N(|gp) > 1, |gp is divisible by some prime vr of o, so g = rp for some r, p G O with |pp an 
associate of vr, by FactlHl Moreover, r and p are reduced, because q is. Now 

(15) qO C Or\qOq-^ C O n rOr^^ 

where the first inclusion is straightforward and the second is a consequence of Lemma |SJ Also 
pOp~^ 7^ O, since p is reduced but not a unit, and hence qOq~^ ^ rOr^^. But O, qOq^^ 
and rOr~^ are maximal orders, so by Lemma |H1 of the Appendix O n qOq~^ 7^ O H rOr^^ 
whence [O n rOr"^ : O n q'C'g""'^]^- / o. On projecting into HI_, (fTK]) gives 

Im(^C') C Im(C' n gOg"^) C Im(C' n rOr"^) = Im(rC'), 

where the equation on the right is an application of the induction hypothesis since r is reduced 
and 

N(|rp) = N(|gp)/N(|p|2) < N(|g|2). 
By Lemma 131 and the multiplicativity of X-indices, we have 

^ ^ ^ ^^^^ [Im(C') : Im(rC')]^ |r|2 

a prime ideal of o, but Fact [7| implies 

[Im(C'nrC'r-^) : Im(C'ngC'g"^)];^ = [O r\rOr-^ : O r\ qOq-^]^ / O. 

Being a divisor of the prime ideal vro, this index must therefore be no itself and, again 
by the multiplicativity of K-indices, [Im(C' n qOq"^) : Im(gC')]^ = o. Hence Im(gO) = 
lm(0 n qOq~^), completing the induction. □ 

Summing up: 

Proposition 4. If K and M{K) both have class number 1, ifO is a maximal order of M{K), 
and if q,r £ O are O-reduced, the following inclusions are equivalent: 

(i) Im(gC') C Im(rC'); 

(ii) qO C rO; 

(iii) OnqOq-^ C ODrOr-^; 

(iv) Im(C' n qOq-^) C Im(C' n rOr~^). 

Proof. Lemma [S] gives (i) =^ (ii) =^=- (iii), the claim (iii) =^ (iv) follows by projection into H^, 
and (iv) <^=^ (i) by Lemma El □ 

As a matter of independent interest, we list the following consequences of Lemma EJ 

Proposition 5. Let M{K) be the standard quaternion algebra over the real algebraic number 
field K , and suppose that both K and M.{K) have class number 1. Let O be a maximal order 
of Il{K) and q £ O an O-reduced element. Then, one has the following relations: 
(i) lm{OnqOq~^) = lm{qO) = lm{Oq); 
(h) O n qOq~^ = O + qO + Oq = O + qO = O + Oq; 

(iii) [0:0n qOq-^] = [lm(0) : lm(0 n qOq-^)] = N(|g|2) ; 

(iv) [Im(C') : lm{0) H lm{qOq-^)] = N(|g|2); 
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(v) [O n qOq-^ : qO] = [o + qO : qO] = [o + Oq : qO] = N(|g|2) . 

Proof. Part (i) follows from Lemma El together with the observation that Im(gx) = lm{—xq) 
when X E O. Part (ii) follows from the inclusions 

O + qO C O + qO + Oq C O H qOq'^ 

and the fact that, by Propositions 121 and Lemma El 

[O n qOq-^ : O + qO]p: = [o : o]^ [Im(C' n qOq'^) : ImiqO)]^ = O. 

For Part (iii), we note that 

[0:OnqOq-% = [o : o]^ [lm(0) : lm(0 n gOg-^)]^, 

by Proposition |2 and use Lemmas IHl and |21 and Fact |S1 Part (iv) follows from Part (iii) 
and Lemma El Finally, Part (v) is a consequence of Parts (ii) and (iii). Lemma the 
multiplicativity of i^-indices and Fact |H1 □ 

6. General results 

We are finally in a position to state our main result on the connection between coincidence 
site modules and one-sided ideals. 

Theorem 1. Let M{K) be the standard quaternion algebra over the real algebraic number 
field K, and suppose that both K and M.[K) have class number 1. Let O be a maximal 
order of M{K). Then, the CSMs of F := Im(C') are in one-to-one correspondence with the 
one-sided ideals qO that are generated by the O -reduced elements q £ O. The corresponding 
coincidence indices are given by N(|gp), so that these right ideals are to be counted according 
to the square root of their ideal index from Proposition (ii) . 

Proof. Parts (ii) and (iv) of Proposition jll provide the connection between the CSMs and the 
right ideals generated by O-reduced elements of O, and Proposition (iii) gives the indices 
of the CSMs. □ 

This implies that we may count the CSMs by counting one-sided ideals in O, which can 
be done by counting all ideals of a given index in O and subtracting off the two-sided ones. 
While the former is achieved by the zeta function Coi^)^ latter is done by the zeta function 
of the base field, via C/^(4s), up to finitely many correction factors from ramified primes. The 
factor 4 in the argument of (j^ stems from the index formula of Lemma ^ Let us introduce 
a Dirichlet series generating function for one-sided ideals that are generated by O-reduced 
elements as Co^('5)j and denote the Dirichlet series generating function for the two-sided ideals 
of O by ^0 0(5), see below for examples. 

Lemma 7. If O is a maximal order in the quaternion algebra M{K) of class number 1, one 
has the identity Coi^) ~ Co'^i^) Co o(^)- 

Proof. Due to the assumption on the class number, each right ideal is principal, and thus of 
the form qO for some q £ O. If g has any divisor that generates a two-sided ideal, we can pull 
that factor out to the right to obtain a unique factorization of qO into a reduced and a two- 
sided ideal. Consequently, sorted by their index, we can either count all right ideals together 
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(giving the term on the left), or reduced right ideals and their possible right multiplication 
by two-sided ideals separately (giving the product of Dirichlet series on the right). □ 

Now, in view of Theorem ^ counting CSMs is the same as counting reduced right ideals, 
but with respect to the square root of the ideal index. This results in 

Theorem 2. Under the assumptions of Theorem^ the Dirichlet series generating function 
<^(s) for the number f{m) of CSMs of index m of the o-module F = Im(C') C M'^ is given by 



m>l 



where Ck(^) '^^ Dedekind zeta function of the field K and E{s) is either 1 or an additional 
analytic factor {with finitely many terms) that takes care of the extra contributions from 
ramified primes. 

In particular, the arithmetic function f{m) is multiplicative, and the elementary coinci- 
dence spectrum of F is U{OC{r)) = i7(S0C(r)) = {m G N | /(m) ^ 0}. 

Proof. The first claim follows from Theorem ^ together with Lemma [3 Since both zeta 
functions involved, C,q{s) and Qq q{s), possess an Euler product decomposition, this applies to 
<?(s) as well, with their detailed structure following from and [20]. As a consequence, f{m) 
must be a multiplicative arithmetic function. The claim about the spectrum is obvious. □ 

Corollary 2. The asymptotic growth of the multiplicative arithmetic function f{m) of The- 
orem is given through its summatory function via 

2.2 

F{x) := /(w) ~ f?"^) as a; —> oo, with g = veSg^2^i^) ■ 

m<x 

In particular, the average size of f{m) is gm. 

Proof. Consider the function ^(s) from Theorem |21 By the general structure of the zeta 
functions involved and the fact that E{s) is an entire function, ^(s) is analytic in the open 
half-plane {s = a + it \ a > 2}, and has a first order pole at s = 2, but no other singularity 
on the line {a = 2}. In fact, ^{s) is analytic everywhere on this line except at s = 2. 

Delange's theorem, compare [4 , Appendix] for a formulation tailored to this situation, now 
gives F{x) ~ gx'^/2, as x — > oo, with g as claimed. □ 

Remark 3. The residue in Corollary [21 can be calculated as 
(16) res,=2^(s) = :^|K^ res,=i C;.(s). 

and reSj,^]^ Ck(^) expressed in terms of values of L-series. 

7. Specific results 

There are three explicit cases of particular interest for applications in crystallography, 
where the ring of integers o is a principal ideal domain (compare Table 1 on p. 422 of [Hj) and 
the quaternion algebra M.{K) has class number 1 (see (22 and the first table on p. 156 of [23] 
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K 


Q 


Q(v^) 


Q(V2) 


O 


J 


I 


K 


o 


z 


Z[r] 





Table 1. Data for three quaternion algebras M.{K) of class number 1, with 
maximal order O. Here, r = (\/5 + l)/2 is the golden ratio, and the ring of 
integers O is a PID in all three cases (see text for details). 



for further examples). The basic information is summarized in Tabled The corresponding 
zeta functions are the Dirichlet series generating functions for the right ideals, 

Ot^SIcO ' n>l 

where 21 runs through the non-zero right ideals of O and /^(m) is the number of such ideals 
of index m. Due to the class number being 1, all ideals are principal, so we have 21 = aO 
for some a S O', and then [O : aO] = N(|a|^) by Proposition |21 (ii). The zeta function for 
the left ideals is the same, as left and right ideals are in a natural one-to-one relation that 
preserves the norm. Similarly, Co o(^) defined as the Dirichlet series that runs over all 
non-zero two-sided ideals of O, so that Co oi'^) ~ Sn>i 

In all three cases of Tabled these zeta functions can be expressed in terms of the Dedekind 
zeta functions of K as follows, compare PH] : 



(17) Co(s) = CKi'^s)CKi2s-l) 



2i-2«), ifi^ = Q, 
otherwise. 



where the extra factor for K = <Q results from the rational prime 2 being ramified in this 
case. The Dedekind zeta functions entering are given by 



c^is) = cis) = riiz^ 



p=±l (5) ^ ^ ^ p=±2 (5) ^ 

'^Q(v^)('^) = 1 _ n d _ r,~s\2 n 1 _ r)-2s ' 

p=±l (8) ^ ^ ' p=±3 (8) ^ 

see |1] for further details in this context, including explicit terms and asymptotic properties. 

We shall also need the Dirichlet series generating functions for the two-sided ideals of O. By 
the relation between prime ideals in O and o according to |2()[ Thm. 22.4], it can be expressed 
by the zeta function of the base field possibly up to finitely many correction factors for 
primes that ramify in the extension from K to M.{K). With the examples of Tabled this 
happens only for J, where (1 -|- i)J = J(l + i) = J(l + i)J is the unique two-sided prime ideal 
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over the rational prime 2. We know from Proposition 13] (ii) that (1 + i)J has index 4 in J, to 
be compared with index 16 for the ideal 2JI. In ah other cases, if ^ is the unique two-sided 
prime ideal in O that lies over the prime ideal p G o, compare |2Ul Thm. 22.3], one has 

Consequently, we obtain for our three examples of Table ^ 

(19) Co.o(s) = Ci^(4s) • < ^, . 

II, otherwise. 

With these general preparations, we can now apply the results to the three cases of Tabled 



7.1. The Hurwitz ring J. Here, we consider EI(Q), with J as defined in Eq. ((71). This is 
the unique maximal order that contains £ ~ of ©, with o = Z and [J : £] = 2. It 
is immediate that Im(JI) = /bco so that we can re-derive the results from Section |2j The 
characterization of the J-reduced elements of JJ follows from the detailed analysis in • 

Fact 10. An element q £ ^ is ^-reduced if and only if q is ^-primitive and {ql"^ is odd. □ 

Theorem 121 applied to J, reproduces all results of Section l2j In particular, we find 

Corollary 3. The elementary coincidence spectrum of the lattice Ibcc = Iiii(J) is the set of 
odd integers. This set is a monoid, generated by the odd rational primes. □ 

In view of Fact llfll one can also derive a simple formula for the index of a rotation matrix 
R G SOC(r'bcc)- R = Rq with g G J being J-reduced, one simply has Uy^cdRq) = kP- 
Otherwise, replace q hy q' = g/contj(g), which is J-primitive. If is odd, the quaternion 
q' is already reduced. If not, one has q' = q"{l + i) with {q"]"^ odd. In this case, since 1 + i 
generates a two-sided ideal, defines the same CSL as Rq'. For a general g G J, we thus 
have 

, V . s I , , s i2 111 if Iq'P is odd, 

(20) I^bcciRq) = |g/contj((7)r- <^ 1 

[21 otherwise. 

This formula is well known and has been derived before |131 0] by different methods. 

The dual lattice of /bcc is the so-called face centred cubic lattice, conveniently written as 

Acc = {ii,m,nY \ i + m + n even} = (l-hi)JnH_. 

It is a sublattice of of index 2. Clearly, due to commensur ability, all three lattices share 
the same commensurator subgroup of S0(3,M), namely S0(3,Q). More generally, we have 

Theorem 3. Let R G S0(3,Q). Then, the three cubic lattices in share the same index 
formula, i.e., Z'bcc(-R) = ^p(-R) = ^fcc{R)- If^ particular, the coincidence spectrum for all 
three lattices is the set of odd integers, as given in Eq. and the Dirichlet series generating 
function for all three cubic lattices is ^(,ub('^) of Eqs. © and (j^J. 
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Proof. By (21 Thm. 2.2], two mutually dual lattices share the same index formula, which then 
applies to Ibcc and Ifcc- So, we have Z'bcc(^) = ^fcc(-R) =: S{R)- 
Moreover, we have the inclusion 

2 32 
/fee C Z C Tbcc • 

By [21 Lemma 2.6], this relation implies Up{R) \ 2U{R) and S{R) \ 2Up(R), which tells us 
that the coincidence index for Z'^ could only differ from I^{R) by a factor of 2. 

Prom our previous arguments, we know that S{R) is always odd, so that I^{R) \ ^p{R) 
On the other hand, also the index i7p(i?) is always odd, see 0]. which emerges from an 
independent argument directly on the basis of the rotation matrix R. Consequently, we must 
have Up{R) \ E{R), which gives Sp{R) = E{R). □ 

Remark 4. Observing that the lattice is isomorphic with J n EI_, the index relation 
Ep{R) = S{R) in the previous proof can also be derived with an independent argument 
based upon Fact El and an index formula for intersections with E[_, see fO] for more. 

Also of interest is the analysis of multiple intersections. Unfortunately, the relation to 
the arithmetic of J becomes considerably more involved. First explicit results are available 
|2ni I2ZI) including the observation that the total coincidence spectrum is identical to the 
elementary one. For a given index, new CSLs appear, but all possibilities are exhausted with 
triple intersections. 

7.2. The icosian ring I. Here, we consider H (Q ( ^5 )) . In this case, there are two maximal 
orders that contain £ of with o = Z[r]. One choice is 

I := ((1,0,0,0), (0,1, 0,0), 1,1,1), i(l-T,r,0,l)>^j^j, 

where r = (\/5 + 1)/2 is the golden ratio. The other choice is F, where ' denotes the algebraic 
conjugation of Q(\/5), defined by ^/b 1— > — \/5- Both have been called the icosian ring., see 
|1()| 1^ ^2 for details on the connection to the root system of type H^. Our choice here 
matches that of jlOl lllj , while the roles of I and F are interchanged in comparison to [2115 . 
It is not difficult to check that £ of © is an index 16 submodule both of I and F. Moreover, 
one has the relation 

2(1 + F) c C c inF, 

where I Pi F = JI[t] := J + rj provides a connection to the Hurwitz ring. Let us mention that, 
in the icosian case, 21 is the unique prime ideal over 2 G Z[r], and that (1 + i)I = F(l + i), so 
that (1 + i)I is a one-sided, but not a two-sided ideal. This also shows that I and F are related 
by an inner automorphism, and hence are of the same type, in line with Fact O Unlike in the 
case of the Hurwitz ring, the rational prime 2, which is also a prime of Z[t], does not ramify. 
More generally, no prime of Z[r] ramifies in the extension step to 1. Consequently, one has 



Fact 11. An element q ^1 is I-reduced if and only if it is I-primitive. 



□ 
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With this definition, the imaginary part satisfies Im(I) = ^^AB, where 

A^B = {Yl^=i ^i^i I ^ with t'^Ui + + a3 = (mod 2)} 

= ((2,0,0)*, (l,l,l)*,(r,0,l)*>^j^, C 

is the standard body centred icosahedral module of quasi-crystahography, see [2j and refer- 
ences therein. Furthermore, I H ]HI_ ~ Im(I H ]HI_) = ^A^F) where 

M.Y = {Yl\=i ^i^i I '^i ^ ^M; with ttj^ = = T^Q^s (mod 2)| 
= {I]Li ^ -^B I + a2 + «3 = (mod 2)} 
= ((2,0,0)*, (r + l,r,l)*,(0,0,2)*>^j^j C A^B C 

is the standard face centred icosahedral module. Both A^b and A^f are free Z[r]-modules of 
rank 3, and free Z- modules of rank 6. One quickly checks that 

[Im(I) : Im(£)] = [I : £] = 16 , 

in line with Fact[7| Note that Im(I) can be obtained as the projection of the weight lattice Dg 
into a 3-dimensional subspace that is invariant under the symmetry group of the icosahedron. 
This also relates A^b to a lattice in 6-space, which is complemented by the corresponding 
relation between AIf and the root lattice D^, see P] for details. 
An application of Theorem [21 to I gives the generating function 

<f, ( \ - 1 + 5" rr /^l±zllV "TT i + p"^' 

_ ^ 5 6 10 24 20 40 30 30 60 

Remark 5. The generating function ^ico('S) applies to both A^b and Alp; as well as to any 
image of either under a linear similarity. The relation between AIb and AIf can be understood 
from the point of view of mutually dual embeddings as lattices into 6-space mentioned before, 
compare (2]. 

The asymptotic growth constant q of Corollary [21 takes the form 

, = res..2^ico(.) = M^^^Ml) ^ 497039^ 

and gm is the average size of the arithmetic function f{m); details of the calculation follow 
from 0J Appendix]. 

The index of a matrix R G SOC(A^b), written as R = Rq with q is given by 
(22) EiRg) = N(|g/conti(g)|2). 

It is thus of the form N(A; -|- £t) = k"^ + k£ — l"^ with k,i E Z. Since each prime vr £ Z[r] has 
a decomposition tt = pp in I, one sees that all indices of this form indeed occur. 
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Corollary 4. The 'Z[T]-modules Mb and Mp have the SOC-group SO(3,Q(a/5)) and share 
the elementary coincidence spectrum 

i:(SO(3,Q(V5))) = {me-N\m = k'^ + M-f, withkJeZ} 

= {1,4,5,9,11,16,19,20,25,...}. 

This spectrum is a monoid and consists of all natural numbers that have all primes congruent 
to ±2 {mod 5) in their factorizations occuring with an even power. 

Proof. The claim on A^b is clear by Theorem |2 as is the characterization of the natm'al 
numbers representable by the quadratic form k"^ + ki — i'^ , compare |14j . 

The SOC-group of M-p is the same because the two modules are commensurate. Also, the 
index formula is the same because the modules Aip and M.^ are dual to one another via the 
embedding into 6-space (as -Dg and Dq, respectively). Consequently, one can use the result 
from [2| that they share the index formula, thus giving equal spectra. □ 

A related object of interest is Z[t]^ = Im(>C), which is a submodule of Im(I) of index 
16 and hence commensurate. Clearly, OC(lm(£)) = OC(lm(C')), and they also share the 
elementary coincidence spectrum 0, though a given index may arise from different sets of 
rotations. Consequently, one gets different Dirichlet series, too, see Sec. 5.3] for details. 

7.3. The octahedral (or cubian) ring K. Here, we consider El(Q(\/2)), with maximal 
order 

117 /l l±i i±j 

^ \-^' V2' V2 ' 2 /Z[V2] ' 

see |1] for a more detailed description. Note that IK contains of ©, with o = Z[-v/2], as 
an index 16 submodule and that K is the only maximal order to contain C In this case, the 
imaginary part is 

Im{K) = -i=((l,0,0)*, (0,1,0)*, ^(1,1,1)*>^[^,, 

so that [Im(K) : Im(£)] = [K : £] = 16, in line with Fact [Tj In the case of K, the rational 
prime 2 ramifies in the extension from Z to Z[\/2], but \/2 does not ramify in the following 
extension to K. As no other prime ramifies, one finds 

Fact 12. An element g E IK is M.-reduced if and only if it is 'K-primitive. □ 
This time, an application of Theorem |21 gives 

^ 1 + 2^' TT (^+P"\' TT 

^oct[s) - Y^^s ii i ) ii 1 _ „2(l-s) 

(23) P=±l (8) ^ ^ ^ p=±3 (8) ^ 

3 6 16 12 10 48 24 36 30 

+ ^ + ^ + ^ + ^ + ^ + Yi^ + Ye^ + Yt^ + Ys^ + • • • 

Here, the asymptotic growth constant g of Corollary [21 reads 

720\/21og(l + \/2) 
Q = res,=2^oct(5) = ' ^ 0.837559, 

with gm again being the average size of f{m). 
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The index of a matrix R = Rg with g G K is given by 
(24) U{R,) = N(|g/contK(<z)H. 

It is thus of the form N(A; + £V2) = k'^ - 2^^ ^ith k,eeZ. Since each prime vr G Z[^/2] has 
a decomposition vr = in K, one sees that all indices of this form indeed occur. This gives 

Corollary 5. The elementary coincidence spectrum of the Ij[^/2]-module Im(]K) is 
Z'(SO(3,Q(\/2))) = {m£N\m = k'^-2f, withk,leZ} 

= {1, 2, 4, 7, 8, 9, 14, 16, 17, 18, 23, 25, . . .} . 

The spectrum is a monoid and consists of all natural numbers that have all primes congruent 
to ±3 [mod 8) in their factorizations occuring with an even power. □ 

A related module in 3-space is Z[\/2]^, realized as Im(£) with the order C of The 
treatment of this module parallels the situation already met with the icosian ring, see above. 

8. Extensions and outlook 

The picture for the coincidence site modules in 3-space that emerge from single intersections 
seems to be rather clear now. Still, there are several important questions that have not been 
addressed here at all. First of all, having counted the CSMs, one would like to have a finer 
classification into Bravais types. Some results in this direction, for the cubic lattices, are 
presented in j25j . Then, there is no compelling reason to stop at single intersections, and 
recent developments make an extension to multiple coincidence site modules desirable. While 
this is rather well understood in the planar case 0, by algebraic means based upon the 
arithmetic of cyclotomic fields, only first steps exist in 3-space, see PHI I2Z| for some results. 

As we briefiy mentioned in passing, several of our above findings can be generalized beyond 
the class number 1 situation, and we hope to report on that soon ^6^ . Quaternions are certainly 
also helpful in 4-space, and further progress in this direction is in sight [7j, at least for simple 
coincidences. More complicated, however, seems the situation in higher dimensions, even for 
the class of root lattices, and that might be a good problem to solve. 

Appendix 

In the proof of Lemma 101 we needed a global version of a uniqueness result on intersecting 
maximal orders in quaternion division algebras, which is to some extent implicit in |2U1 123j . 
but which we could not find in explicit form and which does not seem to have an analogue in 
the more general context of central simple algebras j2()|l24j. For completeness, we sketch here 
a localization argument for this [221 124j that was communicated to us by U. Staemmler. For 
comparison with other work, we recall that the intersection of two (not necessarily distinct) 
maximal orders is called an Eichler order. 

Lemma 8. |22j Let Wi{K) he the quaternion division algebra over the real algebraic number 
field K, and let M, N and O be three maximal orders of IH(i^) that satisfy OCiM = OCiN. 
Then M = N. 
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Proof. We write H = M{K) for this proof. One has M = if and only if ah local completions 
of M and N at prime ideals p oi K are equal, i.e., Mp = Np, see |2S1 Prop. III. 5.1] and 
Sec. 3] for details in the present setting. Note that E := O D M = O D N is an Eichler order. 
Local completion preserves intersections and the property of being a maximal order, see |231 
p. 84]. In particular, Ep = Op D Mp = OpCi Np remains an Eichler order in Hp. 

Either Hp is a skew field (when p divides the discriminant of H over K) or else it is 
isomorphic to the matrix ring Mat(2, i^p) by Wedderburn's theorem, see |17| Sec. 52]. In the 
latter case, Mp = Np by j231 Lemme II. 2. 4], which states that Eichler orders uniquely specify 
the two intersecting maximal orders here. In the former case, Mp = Np because Hp has a 
unique maximal order by |23| Lemme II. 1.5]. So, Mp = Np for all prime ideals p and hence 
M = N. □ 

Note that this proof actually shows a slightly stronger statement: if L, M, N, O are maximal 
orders with L D M = N f] O, one has {L, M} = {N, O}. But the lemma as stated is enough 
for our purposes. 
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